Controlled deflection of cold atomic clouds and of Bose-Einstein condensates 



0^ 
O 
O 
(N 

Oh' 
OO 



Oh; 



> 

m 

0^ 
O 

O 



X 



N. Gaaloul,^'^ A. Jaouadi,^'-^ L. Pruvost,'^ M. Telmini,^ and E. Charron^ 

' Laboratoire de Spectroscopie Atomique, Moleculaire et Applications, 
Departement de Physique, Faculte des Sciences de Tunis, 
Umversite Tunis El Manar, 2092 Tunis, Tunisia. 
^ Laboratoire de Photophysique Moleculaire du CNRS, 
Universite Paris-Sud, Bdtiment 210, 91405 Orsay Cedex, France. 
^ Laboratoire Aime Cotton du CNRS, Universite Paris-Sud, Bdtiment 505, 91405 Orsay Cedex, France. 

We present a detailed, realistic proposal and analysis of the implementation of a cold atom 
deflector using time-dependent far off-resonance optical guides. An analytical model and numerical 
simulations are used to illustrate its characteristics when applied to both non-degenerate atomic 
ensembles and to Bose-Einstein condensates. Using for all relevant parameters values that are 
achieved with present technology, we show that it is possible to deflect almost entirely an ensemble 
of ^'^Rb atoms falling in the gravity fleld. We discuss the limits of this proposal, and illustrate its 
robustness against non-adiabatic transitions. 

PACS numbers: 37.10.Gh 



I. INTRODUCTION 

Optical and magnetic fields are extremely efficient 
tools used for the controlled manipulation of large en- 
sembles of cold atoms [l|, In the past fifteen years, 
cold matter waves have shown great possibilities in the 
context of linear atom optics, when phase-space densi- 
ties are sufficiently low that the effect of collisions can 
be neglected. Dipole and radiation-pressure forces have 
for instance allowed the achievement of various optical 
manipulations such as atomic focusing, diffraction or in- 
terference [1, 13] ■ 

Many efforts have been recently devoted to the exper- 
imental implementation of atomic beam splitters with 
magnetic ll S 0, 111 or optical H, [13, IHl potentials. 
These different experimental investigations were accom- 
panied by various theoretical studies [H, [H, [13, [H, [l^ . 
These devices are obviously of clear interest for atom 
intcrferometry experiments. After the advent of Bose- 
Einstein condensation (BEG) in 1995 [13, El, different 
setups were designed in order to split and recombine a 
BEG [l^, H^, [2l|. In this case, the experimental imple- 
mentation is even more difficult since inter-atomic inter- 
actions due to high atomic densities in the wave-guides 
can sometimes not only induce the fragmentation of the 
BEG [22j2^] , but also affect the overall coherence of the 
system |23|. 

In a recent paper we have derived a semi-classical mo- 
del for the description of the splitting dynamics of a cold 
atomic cloud in such a device (Tsj . This setup involves 
two crossing far off-resonant dipole guides [see Fig.[TJa)], 
and we have shown that a simple variation of the laser 
beam intensities allows to control the splitting ratio in the 
two guides. In the present paper, we first show that if 
the vertical guide is switched off when the atomic cloud 
reaches the crossing point, this device becomes an effi- 
cient coherent atom deflector. We then extend this study 
to the quantum degenerate regime, in order to demon- 
strate the efficiency of this deflection setup with Bose- 




FIG. 1: (a) Schematic representation of the proposed optical 
deflector for cold atoms. The right inset is a magnification 
of the crossing region. The vertical position of the crossing 
point is z = —h, and the total transverse width of the oblique 
guide is equal to 2ii . (b) Timing of the magnetic and optical 
trap (MOT) and of the vertical (V) and oblique (O) guides 
used in this setup, tc = ^/2h/g corresponds to the time at 
which the Rb atoms reach the crossing height z = —h. 



Einstein condensates. 

As illustrated in Fig.[lja), we use a setup involving 
two crossing far off-resonant dipole guides similar to the 
one of Ref. |9| . A large ensemble of ^^Rb atoms is initially 
trapped and cooled around the position z = in Fig.[lja) . 
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This trap is switched off at time t = 0, while a vertical 
far off-resonant laser beam, crossing the cloud close to its 
center, is switched on. A significant portion of the atoms, 
falling due to gravity, is captured and guided in this ver- 
tical wave-guide [Si] . When the center of the guided cloud 
reaches a given height z = — /i, at time t = ic, the vertical 
laser beam is switched off while a second oblique guide is 
switched on. This timing sequence is illustrated schemat- 
ically in Fig.[ljb). The durations of the switching-on and 
-off procedures are supposed to be much shorter than the 
typical time scale of the fall dynamics. In spite of the high 
velocities achieved in this vertical fall, we will show that 
this setup allows for the implementation of an efhcient 
deflector since the atoms can be deviated from their ini- 
tial trajectory with no significant loss. This scheme is 
used both with a thermal cloud of atoms and with an 
atomic condensate after rescaling the whole problem due 
to the difference in size of condensates compared to cold 
atomic clouds. 

The outline of the paper is as follows: in Sec.|TT] we 
discuss the properties of ^^Rb atoms that are relevant 
for our analysis. We also give the values of typical laser 
parameters that realize this atom deflector. We describe 
briefly our semi-classical numerical model in Sec. lIIII In 
Sec. lIVI we give the results of our numerical investigations 
on the performance of this setup with cold atomic clouds 
(T ~ 10 /iK). We show that a high efficiency 90%) 
can be achieved with large deflection angles. We also 
discuss the adiabaticity of the deflection process. We 
then present in Sec.|V]a full quantum model designed to 
treat the dynamics of a BEC falling in the gravity field in 
the presence of these time-dependent guiding potentials. 
We then present the results of the numerical simulations 
with BECs, demonstrating the efficiency of the proposed 
setup in the quantum degenerate domain. Our conclu- 
sions are finally summarized in the last section. 



sian-like form 

if|.K4: /o(x) = _^sm (^-^j , 

if 1x1 >4 : h{x) = 0, 

where the size of the vertical guide is simply related 
to the laser waist vjq by the relation 

4 = woV21n2- 1.18 Wo- (3) 

This sinus-squared shape, which is often used in time- 
dependent calculations [23 , is very close to the ideal 
Gaussian intensity distribution, except for the absence 
of the extended wings of the true Gaussian shape which 
lengthen the calculations without noticeable contribution 
to the physical processes. With this sinus-squared con- 
vention, the guiding region (|x| < ^0) is also well defined. 
The trapping potentials associated with the vertical and 
oblique laser guides are thus expressed as 

U,,{:x) = -t/osin^ (f for \x\ ^ £0 (4a) 

Ui{x,z) = -t/isin^ ^-J^^ for < li (4b) 

where x' denotes the rotated coordinate x' = x cos 7 -t- 
{z + h) sin 7. 

Typical laser powers Po ~ 5 — 30 W for a Nd:YAG laser 
operating at 1064 nm with laser waists of about 100 — 
300 /.tm yield potential depths of about 5 — 250 ^K. With 
these laser parameters, the ^^Rb transition to consider 
is the Di : 5^Si/2 — *'5^Pi/2, with a decay rate V/2tt ~ 
5.75 MHz, a saturation intensity Is — 4.5mW/cm^ and a 
detuning 5/271 ~ — 95.4TIIz. With these conditions, the 
Raylcigh range zr — ttwq/X is about 3 cm, thus allowing 
us to neglect the divergence of the beam on a length up 
to about 1 cm. 



II. GUIDING POTENTIALS 



III. SEMI-CLASSICAL MODEL FOR COLD 
ATOMS 



During the guiding process and in the case of a large 
detuning, the atoms are subjected to a dipole force in- 
duced by the dipole potential 



_ hr I{r)/L 



(1) 



where d — uji^^loq denotes the detuning between the laser 
frequency lol and the atomic transition frequency loq. Ig 
is the saturation intensitj^ and T the natural linewidth 
of the atomic transition [2y, [23| . 

The atomic dynamics is supposed to take place in 
the {x, z) plane defined by the two guides (see Fig.[T](a)) 
thanks to a strong confinement applied in the y-direction. 
The transverse intensity distribution of the TEMqo verti- 
cal laser beam of power Pq is approximated by the Gaus- 



The guided atomic dynamics can be followed by solving 
numerically the time-dependent Schrodinger equation for 
the atomic translational coordinates, taking into account 
the effect of the gravity field, and choosing realistic val- 
ues for all laser parameters. We adopt a semi-classical 
approach where the z coordinate is described classically, 
following 



t>tc: Zci{t) 



(5a) 



-g[tc + {t-tc) con -A /2, (5b) 



where tc — (2/i/g)^/^ is the time at which the atoms 
reach the crossing point (position z = —h). These equa- 
tions of motion are obtained under the assumption of 
energy conservation for a classical particle which is per- 
fectly defiected, and which therefore follows the paths 
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blazed initially by the vertical beam and later on by the 
oblique guide. The other dimension x is treated at the 
quantum level. This semi-classical approach was com- 
pared to the experimental study [9*1 in Ref. . In this 
approach, the two-dimensional guiding potentials (|4]) can 
be replaced by the one-dimensional time-dependent po- 
tential 

t^t^: U{x,t) = Ua{x) (6a) 
t>tc: U{x,t) = Ui{x,zci{t)), (6b) 

and the quantum dynamics is now summarized in the 
one-dimensional time-dependent Hamiltonian 



(7) 



where m denotes the *^Rb atomic mass. The time-de- 
pendent Schrodinger equation 



i-ri—ipix.t) 



n{x,t) (pix,t) , 



(8) 



is then solved using the numerical spli t operator tech- 
nique of the short-time propagator [2g| , assuming that 
the atom is initially {t = 0) in a well defined eigenstate 
V, of energy E^, of the potential (|4ap created by the verti- 
cal laser beam. In addition, it was shown in Ref. 15| that 
the deflection probability obtained for the initial classi- 
cal conditions z{0) = and i(0) = is very close to 
the probability averaged over the entire atomic cloud. 
We therefore use these initial classical conditions in the 
present study. 

At the end of the propagation, the final wave function 
tp{x,tf) is analyzed spatially, in order to extract the de- 
flection efficiency Ty^. An averaging procedure over the 
set of all possible initial states finally allows to calculate 
the total deflection probability (7713) of the entire atomic 
cloud (see Sec. lIV B] hereafter for details). 



IV. NUMERICAL RESULTS FOR COLD ATOMS 
A. Case of a single initial state 

In this study, the value of the position h of the cross- 
ing point between the two guides is the main parameter 
which controls the efficiency of the deflector. Indeed, for 
large values of h the atoms reach the crossing point with 
a large kinetic energy Ec — mgh, and they will not be 
deflected if this energy exceeds by far the binding energy 
in the oblique guide. 

In order to predict precisely the largest value of the 
height h allowing for atomic deflection, one should com- 
pare the kinetic energy gained by the atoms along the 
direction x' transverse to the oblique guide at the posi- 
tion z = —h — £1/ sin 7 [see Fig.[IJa)] with the binding 
energy Ui^Ey. The energy E^ denotes here the energy of 



the initial vibrational state v. One can effectively expect 
that the deflection will fail if 

mg ( h + sin^ J > Ui- Ey . (9) 
\ sm 7 / 

In this expression, the sin^ 7 factor originates from the 
fact that the transverse direction x' of the deflecting 
beam makes an angle 7 with the fall direction z. The 
validity of this simple prediction is illustrated in Fig.[2l 
which represents the deflection probability 77D as a func- 
tion of h [Fig.E^a)] and of wi = £i/V21n2 [Fig.lljb)], 
all other parameters being flxed. These probabilities are 
calculated numerically for the initial state v = and for 
V — 2094, whose energy is about halfway in the opti- 
cal potential {E^ ~ —Uq/2). In both graphs, the fron- 
tiers defined by the inequality ^ arc indicated by ver- 
tical dashed arrows. By comparison with the "exact" 
value obtained from the solution of the time-dependent 
Schrodinger equation ([5]) , one can notice that these fron- 
tiers correspond to a defiection probability of 50%. This 
energy criterion, which simply compares the atomic ki- 
netic energy with the binding energy in the oblique guide, 
can thus be used safely to predict the efficiency of this 
setup. 

One can also notice in Fig. [21^ a) the different variations 
of rjD with ft, for w = and for v — 2094. The differ- 
ent behavior of these two vibrational levels comes from 
the fact that u = is associated with a well localized 
atomic wavefunction, deeply bound in an almost har- 
monic potential, while v — 2094 is entirely delocalized 
over a large spatial range < ^o/2, since its energy is 
about halfway in the potential. As a consequence, v = 
satisfies fully the conditions imposed by the Ehrenfest 
theorem [2^ and its evolution can be described classi- 
cally, while V = 2094 shows a quantum behavior. For 
w = 0, as soon as the inequality ([9]) is satisfied, the de- 
flection probability falls to zero, in agreement with the 
usual dynamics of a classical particle. On the other hand, 
the stationnary state v — 2094 can be seen as a coher- 
ent superposition of incoming and outgoing wave packets 
characterized by a rather broad kinetic energy distribu- 
tion of width AEc ~ Uo/2. The packet moving in the 
-l-a: direction will be easily captured by the oblique guide, 
while the packet moving in the opposite direction easily 
avoids this wave guide. These two different dynamics are 
not much affected by the exact value of the falling height 
h, and this explains the very slow variation of 77/3 with h 
in Fig.UJa) for v = 2094. 

The variation of iju with wi [see Fig.^b)] is also op- 
posite for V = and v = 2094. The case v = can again 
be interpreted classically : when wi increases, the possi- 
bility is open for the atoms to fall from a higher distance 
d — {h + £1/ sin 7), thus gaining a larger kinetic energy. 
This explains the decrease of 77D with wi for v — 0. This 
variation is just reversed in the case oi v — 2094. Here, 
the initial wave function is characterized by a large typ- 
ical size Ax ~ £0. An efficient deflector can thus only be 
obtained if the size of the oblique wave guide remains of 
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FIG. 2: Deflection probability rjo as a function (a) of the 
falling distance h [see Fig.[T] and (b) of the waist wi of the 
oblique laser beam. The deflection angle is equal to 7 = 
10 deg. These results are for a single initial state : v = 
(solid line with red circles) or « = 2094 (solid line with green 
squares). In both graphs, the dashed blue arrows mark the 
positions at which a deflection efficiency of 50% is expected 
according to inequality Q . The laser parameters have been 
chosen such that Uo — Ui = 30 nK, and wq — 100 fim. This 
corresponds to P/5 ~ 3.1 x 10~''W/GHz. In graph (a) the 
oblique laser waist is wi = 100 /im and in graph (b) the height 
h is equal to 9.03 mm for v = and to 4.18 mm for v = 2094. 



the order of, or is higher than, this typical size £0. Con- 
sequently, for V — 2094, when wi decreases below wq, the 
deflection probability decreases, as seen in Fig.[2llb). 

In addition, it is worth noting that, due to the sin^7 
factor in the inequality it is possible to induce an 
efficient deflection of atoms with relatively large kinetic 
energies using modest laser powers, as long as the angle 
7 remains small. For instance, in the case v = shown 
in Fig.[2Ia) for 7 — 10 deg, an almost perfect deflection 
is obtained for h = 8.5 mm, even though the total kinetic 
energy of the atom reaches then about Ec ^ 900 /iK, i.e. 
30 times the depth of the oblique wave guide. A larger 
deflection angle could be achieved easily and with a very 
high efficiency by simply adding a succession of several 
deflection setups, each one inducing a small deflection of 
about 10 degrees. 

An important issue for the preservation of the coher- 
ence properties of an atomic cloud is the adiabaticity 
of the process. Previous theoretical studies have shown 
that similar beam splitter setups are able to conserve the 
coherence of the system even for a thermal distribution 
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FIG. 3: Atomic probability distributions \(p{x,tf)f as a func- 
tion of the transverse coordinate x at the end of the prop- 
agation, for (a) h =2 mm, (b) h =7 mm and (c) h = 9 mm. 
The laser parameters are identical to the one of Fig. (2] with 
V — and wi = 100 /im. Note that, for the sake of clarity, the 
horizontal axis has been broken in panel (c). The small insets 
in panels (a) and (b) represent the vibrational distributions 
in the oblique guide at the end of the propagation. 



of atoms with an average energy far exceeding the level 
spacing of the transverse confinement ^IS'l . This behavior 
results from the fact that non-adiabatic transitions are 
induced by the time derivative operator d/dt which does 
not couple states of opposite parities, thus preventing 
nearest neighbor transitions [10] ■ In comparison, transi- 
tions to other states presenting the same parity as the 
initial state are also not favored since they involve larger 
energy differences [l^ . 

As shown in Fig. [31 this robustness to non-adiabatic 
transitions is also present in our deflection scheme. This 
figure represents the probability distributions \(p{x,tf)\^ 
calculated 7 mm below the crossing point z — —h for the 
initial state v — 0, with h — 2 mm [Fig.[3Ia)], h = 7 mm 
[Fig-Elb)], and h — 9 mm [Fig.[3fc)]. The vibrational dis- 
tributions obtained in the oblique guide after deflection 
are also shown in the small insets of Fig. [31^ a) and[3I^b). 
Even though the kinetic energy of the atoms exceeds the 
average vibrational spacing in the trap, the initial state 
w = is preserved at 99.1% for h — 2 mm, and at 50.3% 
for h ~ 7 mm. Indeed, in the first case, only w = 2 is 
slightly populated, while the first five even vibrational 
levels are populated in the second case. It is only when 
the falling height h approaches the limit given by the in- 
equality ^ that the population of the initial state v = 
is almost entirely redistributed to higher excited states, 
as seen in the wave function shown Fig.[3Jc). 



5 



1000 



2000 



3000 



4000 



5000 



• f/j = 30 |iK 

□ f/j = 25 |iK 




(a) h = 4 mm 



5000 



FIG. 4: Deflection probability r;_D as a function of the initial 
vibrational level v. (a) The deflection angle is equal to 7 = 
lOdeg, and the solid line with red circles stands for Ui = 
30 /iK while the solid line with green squares is for Ui — 
25 fjK. (b) The depth of the oblique wave guide is equal to 
Ui = 30 /iK, and the solid line with red circles stands for 7 — 
10 deg while the solid line with green squares is for 7 = 13 deg. 
The falling height is ft = 4 mm and the oblique laser waist is 
equal to wi — 100 /im. All other parameters are identical to 
the one of Fig.[3 
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FIG. 5: Total deflection probability (r/o) of an atomic cloud of 
^'^Rb of size ctq — 0.15 mm at temperature To — 10 jjK. The 
laser parameters have been chosen such that Uo = 30 /iK, 
wq = 200 /im, wi = 158 fim, and h = 4 mm (a) or ft = 1 mm 
(b). 



B. Case of an initial vibrational distribution 

Realistically, an atomic cloud of typical size (Tq and 
temperature Tq can be described as a statistical mixture 
of trapped vibrational states, represented by the density 
matrix 

p(cro,To) = ^c„(cro,7o) \v){v\, (10) 

V 

where the coefficients Cy{ao, Tq) are involved functions of 
the cloud parameters ctq and Tq and of the wave guide pa- 
rameters Uq and Wq (see equation (16) in reference [l^ for 
instance). The calculation of the total deflection prob- 
ability of the entire cloud therefore requires to average 
incoherently the deflection probability of each possible 
initial vibrational level v, taking into account the weight 
functions Cy{ao,To). It is also worth noting that typical 
initial vibrational distributions P{v) — |cu(cro,To)p are 
relatively flat when fc^T ^ Uq, except for the lowest en- 
ergy levels which are usually more populated [lH| • In the 
calculation, we include all populated vibrational states. 

Fig.[3] shows the variation of the deflection efficiency 
with the initial vibrational level v for a series of different 
laser parameters. The transverse trapping potential asso- 



ciated with the vertical wave guide supports about 5000 
vibrational states when Uq — 30 fiK and wq — 100 /im. 
One can notice the general tendency of measuring a lower 
deflection probability when v increases, in perfect agree- 
ment with the variation expected from the energy cri- 
terion®. In addition, one can notice that increasing 
Ui [Fig.HKa)] or decreasing 7 [Fig.[Hb)] increases the de- 
flection probability of any initial state. In Fig.lH the 
vertical dashed arrows indicate the limits defined by the 
inequality®, which are again in good agreement with 
the numerical values. One can also notice that the high- 
est levels V ~ 5000 are not deflected. This is due to the 
fact that atoms trapped in these levels, whose energies 
are very close to the threshold, are easily lost during the 
deflection process. 

Fig. [5] represents the averaged deflection probability 
(?7d) as a function of the deflection angle 7 and of the 
potential depth Ui of the oblique laser guide, for a ther- 
mal input state of size ctq = 0.15 mm and temperature 
To = lO^K, with h — 4 mm [Fig.[Ha)] and h = 1mm 
[Fig-EIb)]. Realistic values have been chosen for all laser 
parameters, close to the one used in the experimental 
study 0, and the coefficients c„(cro,?o) of Eq. pUj) were 
calculated following Ref . [l^ . One can notice a rapid de- 
crease of (?7z)) when Ui decreases and when 7 increases. 
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However, an almost complete deflection (93. 8%) is still 
observed in the case h = 1 mm with 7 = 25 deg and 
Ui = 120 /xK, even though the total kinetic energy of 
the atoms reaches then about Ec ~ 100 /iK at the cross- 
ing point, all trapped states being significantly populated 
initially. For 7 = 10 deg, the deflection efficiency reaches 
99.8%. We have also verified that decreasing the temper- 
ature of the initial cloud increases significantly the deflec- 
tion efficiency since it suppresses the population of the 
highest trapped states, for which the deflection process 
is less efficient [see Fig. 3]. It is also worth noting that 
since the deflection process is less efficient for the highest 
trapped levels, it could also be used to selectively sepa- 
rate the lowest energy levels of the trap. Since it behaves 
very well for the lowest trapped states, we expect that 
this setup will prove useful with Bose-Einstein conden- 
sates. We therefore derive in the next section a quantum 
model aimed at the description of the dynamics of a Bose 
gas in such a deflection setup. 



V. DEFLECTION OF BOSE-EINSTEIN 
CONDENSATES 

A. Theoretical model 

From the theoretical point of view, in the case of a 
low density the dynamics of the macroscopic wave func- 
tion ^'(r, t) of a Bose-Einstein condensate can be accu- 
rately described by the mean-field Gross-Pitaevskii equa- 
tion js^, [31I, [S^l ■ In three dimensions and in the presence 
of both a time-dependent external potential V{r, t) and 
the gravity field this equation reads 



- —V^ + Vir,t) + mgz + NUol'fi 
Zm 



(11) 



where Uq = ^-kIt? / m is the scattering amplitude and oq 
the s-wave scattering length. N denotes the condensate 
number and NUq \^{r,t)\^ is the mean field interaction 
energy. The three-dimensional coordinate is denoted by 
r= {x,y,z). 

In the absence of a trapping potential in the z-di- 
rection, the condensate will not only expand but also 
fall around the average classical height Zci{t) — 
Since the de Broglie wavelength of the BEG is no more 
negligible compared to the characteristic distances of the 
problem, a quantum treatment of this direction is neces- 
sary unlike the thermal atoms case discussed in the first 
part of the paper. The simulation of the fall dynamics is 
thus greatly facilitated when done in the moving frame 
R= iX,Y,Z), where 



1 9 

using the unitary transformation 



S(J?, t) — exp 



^'(r,i). 



(12) 



(13) 



Indeed, applying this transformation yields a simplified 
Gross-Pitaevskii equation 



^Vl + V{R,t) + NU^\E.\^ 



(14) 



where the gravitational term mgz has vanished. 

Following the variational approach of Ref. [s^ , we now 
assume (as in Sec. IIIII and IIV[) a strong harmonic con- 
finement in the perpendicular ^-direction, with 

V{R,t) = ]^mijjlY^ + V\\{X,Z,t) , (15) 

where V\\ {X, Z, t) denotes the optical guiding potential 

t^t,: V\\{X,Z,t) = Uo{X) (16a) 
t>tc: V\\{X,Z,t) = UiiX,Z - gf/2) . (16b) 

The confinement along the perpendicular direction Y is 
supposed to be much stronger than along the parallel 
directions X and Z , thus yielding the conditions 



4C/o 



and uj± 3> 



4Ui 



(17) 



The condensate dynamics is now followed using the ap- 
propriate ansatz [H, [s^ 



S(i2,t) = $(X, Z,t)f{Y\n) 



where 



f{Y\n) = 



1 



(18) 



(19) 



This choice amounts to assume a Gaussian shape of the 
wave function in the l"-direction, characterized by a time- 
dependent width VL{X^Z,t). This width varies slowly 
along the parallel directions, thus implying 



(20) 



It has been shown that this choice is well justified not 
only in the limit of weak interatomic couplings but also 
with large condensate numbers 36 , 37]. 

An effective two-dimensional non-linear wave equation 
is then derived using the quantum least action princi- 
ple (HIH for $(X, Z,t) 



^9$ 



2m ■ I 
1 / 
^4 I mr22 



V2^ ^ 



(21) 



This equation describes the condensate dynamics in 
the {X, Z) plane, with an accuracy which goes beyond 
the usual two-dimensional Gross-Pitaevskii equation. It 
takes into account the influence of the dynamics along 
the perpendicular direction on the evolution of <I>(X, Z, t) 
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FIG. 6: (a) and (c) : Atomic density |<l>(x,2,t)p in arbitrary 
units as a function of x and z at times t = (upper graph) 
and t — 5.2 ms (lower graph), (b) and (d) : Width parameter 
^2(2;, z, t) in /im as a function of x and z at times t = (upper 
graph) and t — 5.2ms (lower graph). In each sub-plot, the 
color scale is defined on the right hand side of the graph. The 
guiding potentials are defined by Uo = 2.2 /xK [P/5 ~ 2.3 x 
lO""^ W/GHz], Ui = 8.8/iK [P/S ~ 9.2 x 10"^ W/GHz], and 
Wo — wi — 300 fim. The deflection angle is 7 = 50 deg and the 
condensate number is TV = 5 x 10*. In the j/— direction, the 
trapping frequency LUy is assumed to be 10 times larger than 
ujx- The crossing height is z = —h — — lO/im and the time 
step for the split operator numerical propagation is St = 1 fis. 

with the introduction of the width parameter ^{X, Z, t). 
Note the difference by a factor 1/2 in the last two terms 
of Eq.dUl) when compared to Eq.(25) of Ref. ^sl, due to 
a misprint in Ref. [33| . The least action variational prin- 
ciple also yields the following quartic equation governing 
the evolution of this parameter 

(^,„.i-^)a.-^l*po-|l.o, ,22) 

\2 w^{t) J 2v27r 2m 

where w{t) = wq for t ^ tc and u'(i) = wi for t > tc- 
This last equation was obtained assuming fl(X,Z,t) <JC 
w{t) for all X, Z and t, in agreement with the strong 
confinement in Y. Compared to Eq.(26) of Ref. [s^, the 
additional term 2V\\/w'^{t) is a small correction due to 
the F-dependence of the TEMqo laser intensity profile. 

The time-dependent wave equation (|2ip is solved nu- 
merically u sing the splitting technique of the short-time 
propagator [2^ , while the quartic equation (|22p for fl is 
solved at each time step and for each coordinate grid 
point {X, Z) using an efficient numerical algorithm [39| . 
The initial wave function is obtained using the imaginary 
time relaxation technique (40| , and the three-dimensional 
condensate wave function ^(r, is reconstructed at the 
end of the propagation by inverting the transforma- 
tions (HH) and p^ . 



B. Numerical Results 

A typical BEC dynamics is illustrated in Fig. [51 which 
shows a surface plot of the atomic density \^{x,z,t)\'^ 
and of the width parameter n{x, z, t) at time t = Q [upper 
part, labels (a) and (b)] and during the deflection process, 
at time t = 5.2 ms [lower part, labels (c) and (d)]. In 
this numerical example, the condensate number is N = 
5 X 10'*. The crossing height z — is reached at time 
tc — 1.43 ms, and the deflection angle is fixed at the value 
7 = 50 deg. In the lower graphs the transparent oblique 
line indicates the direction of deflection corresponding to 
this angle. 

One can first notice in Fig.[BI^b) and (d) that even for 
a relatively small condensate number, a Gaussian ansatz 
for f{Y\Vl) can only be used if the width is a free pa- 
rameter which can take different values at different po- 
sitions X and z. This result is in agreement with the 
numerical studies of Salasnich et al [33| . 

The inset (a) of Fig. [6] shows the symmetric ground 
state wave function of the initial condensate prepared in 
a parabolic trapping potential of equal frequencies uj^ = 
LOz- At time t = 5.2ms [insets (c) and (d) of Fig.[S], the 
condensate wave function has expanded during the fall 
dynamics, and has been efhciently deflected along the 
direction 7 = 50 deg by the oblique laser guide. 

In our numerical simulations, we propagate the con- 
densate wave function well after the crossing point be- 
tween the vertical and oblique laser beams has been rea- 
ched, and we obtain the deflection efficiency rj u hy cal- 
culating the condensate number in the oblique trapping 
potential at the end of the propagation. Fig. [7] shows the 
variation of the condensate deflection efficiency with the 
crossing height h [inset (a)], the deflection angle 7 [in- 
set (b)] and the ratio of laser powers Pi/Pq = Ui/Uq 
[inset (c)]. The other parameters are given in the figure 
caption. 

The simple energy criterium ([9]) is still shown to be 
relatively accurate with Bosc-Einstcin condensates since 
in Fig.[7Ua), a deflection efficiency of 50% is expected 
for a falling height h = 114/j,m according to Eq.®, and 
the numerical simulation yields a deflection probability 
of 55 %. The variation of the deflection efficiency rjn with 
/i, 7 and Ui/Uq is also found to be very similar to the 
one obtained with cold atomic clouds. One can finally 
notice that large deflection angles can be reached when 
Ui > Uq, with almost no atom loss for instance when 
7 = 50 deg and Ui — AUq. This high deflection efficiency 
can be explained by the small size of the BEC and the 
weak spread in velocities compared to cold atoms. 



VI. CONCLUSION 

In summary, we have presented a detailed analysis of 
the implementation of an optical deflector for cold atomic 
clouds and for Bose-Einstein condensates. Our analy- 
sis is quite close to the experimental conditions, and is 
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clearly within the reach of current technology. We have 
shown how to create a high performance deflector us- 
ing two crossing laser beams which are switched on and 
off in a synchronized way. We have found that a 10 /iK 
cloud of Rubidium atoms can be deflected by 25 degrees 
with an efficiency of about 94%, and by 10 degrees with 
an efficiency exceeding 99%. A succession of such de- 
flecting setups at this small angle could also be imple- 
mented in order to achieve larger deflection angles with 
high fidelities. We have shown that this device is robust 
against non-adiabatic transitions, an undesirable effect 
which could have led to heating processes. A high degree 
of control can therefore be achieved with such quantum 
systems, opening some possibilities for a range of applica- 
tions. We have also derived an original approach treating 
the dynamics of a Bose-Eintein condensate in the gravity 
field using the quantum least action principle in a mov- 
ing frame. This model was used to demonstrate the high 
efficiency of this defiection setup with quantum degener- 
ate gases since deflection angles of up to about 50 degrees 
can be implemented with no significant atom loss. 



Acknowledgments 



FIG. 7: Condensate deflection efficiency rjD as a function of 
(a) the crossing height h, (b) the defiection angle 7 and (c) 
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laser powers Pi/Pq. In these three graphs, the potential depth 
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In (b) Ui = 8.8 fjK and h — 10 ^m. In (c) h — 10 /im and 
7 — 50 deg (blue line with circles) and 7 = 70 deg (red line 
with squares). The dashed blue arrow in the upper graph 
marks the position at which a deflection efficiency of 50 % is 
expected according to Eq.(|9]). 
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